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Abstract 

The diffusion equation for nonequilibrium interstitial impurity 
atoms taking into account their charge states and drift of all mo- 
bile interstitial species in the built-in electric field and in the field 
of elastic stress was obtained. The obtained generalized equation is 
equivalent to the set of diffusion equations written for the interstitial 
impurity atoms in each individual charge state. Due to a number of 
the characteristic features the generalized equation is more convenient 
for numerical solution than the original system of separate diffusion 
equations. 

On this basis, the macroscopic description of stress-mediated im- 
purity diffusion due to a kick-out mechanism was obtained. It is sup- 
posed that the interstitial impurity atom makes a number of jumps 
before conversion to the substitutional position. At the same time, 
a local equilibrium prevails between substitutionally dissolved impu- 
rity atoms, nonequilibrium self-interstitials, and interstitial impurity 
atoms. Also, the derived equation for impurity diffusion due to the 
kick-out mechanism takes into account all charge states of interstitial 
impurity atoms as well as drift of interstitial species in the electric 
field and in the field of elastic stress. Moreover, this equation exactly 
matches the equation of stress- mediated impurity diffusion due to gen- 
eration, migration, and dissociation of the equilibrium pairs "impurity 
atom ~ self-interstitial" . 

Keywords: Theory of diffusion; diffusion; exchange interactions; 
annealing; semiconductors 
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1 Introduction 



In recent years a mechanism of dopant diffusion in silicon due to forma- 
tion, migration and dissociation of the pairs "impurity atom - vacancy" or 
"impurity atom - self-interstitial" (pair diffusion mechanism) has become 
commonly accepted (see, for example [H El IS])- However, boron diffusion 
in silicon is often considered within the framework of the substitutional- 
interstitial mechanism U El E] when the silicon self-interstitial displaces an 
immobile impurity atom from substitutional to the interstitial position. A 
migrating interstitial impurity atom in turn replaces the host atom becoming 
substitutional again (so-called "kick-out mechanism"). 

It is supposed that kick-out mechanism is also responsible for the diffusion 
of gold in silicon % zinc [3110111111121113111111131121, nitrogen [HI UHl 
HH, magnesium 120112111221, and beryUium [m [I11I2D1I23 in GaAs. As 
follows from the data of investigations |14t 1^ 1^ . beryllium diffusion in 
other compound semiconductors is governed by kick-out mechanism too. 

A number of diffusion equations were proposed [HllllZl UHl HH UHl UH UHl 
HH 1221 1231 1211 123 I26j to describe migration of interstitial impurity atoms. 
The equations used in jH [71 [TUl [Hj do not take into account the electric field 
effect on the migration of interstitial impurity atoms and can be used only 
for modeling of the diffusion of neutral interstitial species. The influence of 
built-in electric field on the drift of charged impurity interstitials is taken into 
consideration in [H [TTl [TH [THl [211 [23 , but diffusion equations were obtained 
only for the case of a single charge state of impurity interstitials. The general 
case of diffusion of impurity interstitial atoms in different charge states was 
considered in [T31 [221 1231 1201 • In these papers the equations for diffusion of 
impurity interstitials and impurity diffusion due to kick-out mechanism were 
obtained. The sum operation of diffusion equations for impurity intersti- 
tials in different charge states, the mass action law for conversions between 
different charge states of diffusing species, and the mass action law for the 
conversion reactions of impurity atoms between substitutional and interstitial 
sites were used in [T31 [221 [23j to obtain the diffusion equation for impurity 
interstitials. In contrast to [131 [221 [23] , in paper [201 the mass action law 
for conversions between different charge states of diffusing species was used 
to obtain the diffusion equation in which only the concentration of neutral 
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impurity interstitials C^^^ is included. This equation has the form 
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where D^^ ( x) is the effective diffusivity of interstitial impurity atoms; x is 
the concentration of charge carriers (electrons or holes) normalized to the 
intrinsic carrier concentration n,; k'^^i^x) is the effective coefficient of absorp- 
tion of interstitial impurity atoms due to kick-out of host atoms from the 
lattice sites; is the effective coefficient describing absorption of in- 

terstitial impurity atoms due to their interaction with vacancies (can be used 
for simulation of Frank- TurnbuU diffusion mechanism [2Z1); k'^ [x) is the ef- 
fective coefficient describing generation of interstitial impurity atoms due to 
replacement of the impurity by self-interstitials from the substitutional posi- 
tion to the interstitial one (Watkins effect PH]); hf^^ are the constants of the 
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local equilibrium for conversions between different charge states of impurity 
interstitials; k is the charge state of interstitial impurity atom A^; z and z^^'^ 
are respectively the charges of impurity atom in substitutional and intersti- 
tial positions in terms of the elementary charge (z = + 1 in the case of doping 
by donors and z = -1 for an acceptor impurity); C^^ and C-^ are the actual 
and equilibrium concentrations of self-interstitials in the neutral charge state; 
C^^ and C^^ are the actual and equilibrium concentrations of vacancies in 
the neutral charge state; S^'"^ and G^'^'^ are respectively the effective rates 
of absorption and generation of impurity interstitials during the formation 
(annealing) of the clusters of impurity atoms; 5*"^^^ and G"^^^ are respectively 
the effective rates of absorption and generation of impurity interstitial atoms 
due to evolution of radiation defects. As can be seen from the expressions 
fl3l4l5l6p and ()8l9p . the effective coefficients and effective functions are ob- 
tained by sum operation of analogous quantities for each individual charge 
state k. 

Eq.(^ has the following characteristic features, i) In comparison with the 
equations of [131 1221 122] , the obtained equation can describe long-range mi- 
gration of nonequilibrium impurity interstitials because the mass action law 
for substitutional impurity, self-interstitials, and impurity interstitial atoms 
was not used, ii) This equation describes diffusion of all interstitial impurity 
atoms with different charge states as a whole, although only the concen- 
tration of the neutral impurity interstitials C^^^ must be obtained to solve 
the equation, iii) The obtained equation takes into account the drift of all 
charged species due to the built-in electric field. At the same time, there is 
no explicit term describing the drift. 

It is to be noted that Eq.(^ is very convenient for numerical solution 
owing to the features ii) and iii). However, similar to the equations of 
jl3L 1221 122] 1 Eq-© ignores the drift of impurity interstitials in the field of 
elastic stresses. But the available experimental data show that elastic stresses 
can significantly influence the impurity diffusion [.29,. 30, 31^. including the 
case of Be diffusion in p+ GaAs [29 . Thus, the macroscopic description 
of stress-mediated diffusion of interstitial impurity atoms requires further 
investigations. 

Yoshida et. al. [32] have investigated qualitatively the changes in the 
lattice configurations of impurity atom and self-interstitial in silicon for kick- 
out and pair diffusion mechanisms and concluded that there is no essential 
difference between the macroscopic descriptions of both mechanisms of im- 
purity migration. However, a rigorous proof of this statement has not been 
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given. Moreover, in |32I the charge states of impurity atoms and point de- 
fects were not considered, and the drift of charged particles in the built-in 
electric field was not taken into account. It was mentioned in 01221 ESI El] 
that the pair diffusion model is mathematically equivalent to the kick-out 
mechanism, but only in [22] the equation for impurity diffusion obtained for 
kick-out mechanism was compared with the diffusion equation for the pair 
diffusion mechanism. In all these papers the case of stress-mediated diffusion 
was not considered. Besides, it was supposed in [221 that a local equilibrium 
between substitutionally dissolved impurity atoms, self-interstitials, and im- 
purity interstitial atoms prevails. This condition is easily implemented if 
the generated impurity interstitial atom makes only one jump before con- 
version to the substitutional position. But in this case the jump direction 
is dependent on the gradient of self-interstitial concentration. It is supposed 
for the pair diffusion mechanism that the pairs make a number of jumps 
before their dissociation. The plain fact is that the direction of the pair mi- 
gration is independent of the gradient of self-interstitial concentration. The 
diffusion equations for both mechanisms will be different if the impurity in- 
terstitial atom makes only one jump while the distribution of self-interstitials 
is not uniform. Therefore, it seems reasonable to verify the identity of the 
macroscopic descriptions for the pair and kick-out diffusion mechanisms if 
the stress-mediated migration occurs. 

The purpose of this study is to obtain the equation of the diffusion of 
interstitial impurity atoms taking into account all their charge states and 
drift of all interstitial species in the built-in electric field as well as in the field 
of elastic stresses. On this basis, the identity of the macroscopic descriptions 
of both diffusion mechanisms for stress-mediated diffusion can be confirmed. 



2 Diffusion equation for nonequilibrium in- 
terstitial impurity atoms 

Taking into consideration the charge states of impurity atoms and self- 
interstitials, the reaction for kick-out mechanism can be written in the form 
similar to that used in 01 122] 

As + l' + m,e- ^ A^, (10) 
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where As and are the impurity atoms in substitutional and interstitial 
positions, respectively; F is the self-interstitial in the charge state q; nie is 
the number of electrons participating in the reaction. 

It is important to note that the charge conservation law 

z + z'i-m^ = z^ (11) 

is valid for reaction ()lUp. Here z'^ is the charge of a self-interstitial. 

Let us suppose that an interstitial impurity atom makes a number of 
jumps before conversion to the substitutional position, i.e. migration of in- 
terstitial impurity species in the absence of electric and stress fields is chaotic. 
Then, the influence of cross-effects is negligible, and the diffusion equa- 
tion of the interstitial impurity atoms with the charge state k can be written 
in the following form: 



QfjAIk ^ Q /^^^^gQAIk\ ^ Q /^^AIkjjAIk(jAIk q^^ 
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where C^^'" and /^^^'^ are the concentration and diffusivity of interstitial 
impurity atoms in the charge state k, respectively; JJ^^'' is the potential 
energy of the interstitial impurity atom in the charge state k in the field of 
elastic stress; 5*"^^^ and G^^^ are respectively the absorption and generation 
rates of the interstitial impurity atoms in the charge state k per unit volume 
of semiconductor. 

The terms in brackets in the right-hand side of Eq. (jl2p . taken with the 
minus sign, describe the fluxes of impurity interstitials due to the concentra- 
tion gradient, drift in the built-in electric field, and drift under the influence 
of elastic stress, respectively. 

Let us assume that interstitial impurity atoms have different charge states 
including the neutral one. In this system of equations (|12j) written for 

all charge states k should be used for description of the impurity interstitial 
diffusion. 

The conversions of interstitial impurity atoms between different charge 
states are described by the following reactions: 

- z'^'^e- ^ k'^'K (13) 
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Due to high mobihty of electrons (holes) there is a local thermodynamic 
equilibrium between the interstitials in different charge states and electrons 
(holes). Then, for reaction (fTSj) the relation 
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can be obtained from the mass action law. Here h^^^ is the constant of a 
local thermodynamic equilibrium for this reaction. 

Substituting the concentration of impurity interstitials C^^^ from the 
relation ()14|) in Eq. ()12|) . we transform this equation to the following form: 
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local thermodynamic equilibrium in the heavily doped region h^^^ from their 
intrinsic values in the undoped semiconductor hf^^. 

Let us assume that the functions a'^^'' depend weakly on the charge state 
of impurity interstitial, i.e. a"^^^ ~ a"^^ holds. Then, summing the equations 
fll5|) written for different charge states /c, one can obtain the generalized 
diffusion equation for nonequilibrium interstitial impurity atoms 
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In Eq. lfT^ quantities ( x) and v^^ (x, x) a^^e the effective diffusivity of 
interstitial impurity atoms and effective drift velocity of impurity interstitials 
in the field of elastic stress, respectively. Note that the expressions (j21) and 
Q are valid for the quantities ^"^^ ( x) and D^^ ( x) ■ 

Similar to Eq.((TI), the obtained diffusion equation (fTBj) is characterized by 
the following features, i) This equation can describe the long-range migration 
of nonequilibrium impurity interstitials. ii) This equation describes diffusion 
of all the interstitial impurity atoms with different charge states as a whole, 
that is a set of equations (fT^ is replaced by one generalized diffusion equation 
fll6|) . iii) The only variable that should be estimated in this equation is the 
concentration of neutral impurity interstitials. iv) The obtained equation 
takes into account the drift of charged species due to the built-in electric field. 
However, there is no explicit term describing the drift due to the electric field, 
v) Although the effective coefficients of Eq. p6p represent nonlinear functions 
of X ill comparison with the constant transport coefficients of equations ()12|). 
these functions are smooth and monotone, vi) Eq. ()16|) describes in general 
diffusion of the nonequilibrium interstitial impurity atoms regardless of the 
mechanism of interstitial generation. 

Note that due to the above-mentioned features, Eq.lfT^ is more conve- 
nient for numerical solution than a set of equations (fT^. 

3 Diffusion equation for interstitial impurity 
atoms under equilibrium conditions 

Now, on the basis of Eq. (fTB|) . the interstitial diffusion due to a kick-out mech- 
anism can be investigated. Within the framework of the kick-out mechanism 
the impurity interstitial is generated as a result of interaction between the 
substitutionally dissolved impurity atom and self-interstitial via the reaction 
fllO|) or due to generation of the interstitial impurity atom in the neutral 
charge state via the reaction 

As + I^ +m^^^e- ^A]^, (18) 

where Af and are the interstitial impurity atom and self-interstitial in the 
neutral charge states, respectively; m^^^ is the number of electrons partici- 
pating in the reaction ()18|). where from the charge conservation law it follows 
that m^^^ = z. 
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Usually, it is assumed that there is a local thermodynamic equilibrium 
between substitutionally dissolved impurity atoms, self-interstitials, and im- 
purity atoms in the interstitial position. Then, for reaction (fTHj) the mass 
action law can be written in the form 
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where C is the concentration of substitutionally dissolved impurity atoms; 
C^^ is the concentration of self-interstitials in the neutral charge state 
H^^^ is the constant of a local thermodynamic equilibrium for reaction ()18|1 . 
With m^^^ = z relation (fT^ may be as follows: 

C^ix = H^i>^^C^>^C. (20) 

Substitution of (pUjl into (fTB|) gives the generalized diffusion equation of 
interstitial impurity atoms in the following form: 
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Eq. is not so general as Eq.(fTH|). because it gives no description for 
the long-range migration of the nonequilibrium impurity interstitials. On the 
other hand, Eq. (PT|) is valid for the very important case of the local equi- 
librium between substitutionally dissolved impurity atoms, self-interstitials, 
and impurity atoms in the interstitial position. The concentration of intersti- 
tial impurity atoms is not included in the explicit form in Eq. (j21|) . although 
diffusion occurs due to migration of interstitial impurity species in different 
charge states. 



4 Equation of impurity diffusion due to kick- 
out mechanism 

Based on Eq. ()21|) . we can obtain the equation of impurity diffusion due to 
a kick-out mechanism. It is commonly accepted that the substitutionally 
dissolved impurity atoms are immobile. Then, the conservation law for sub- 
stitutional impurity atoms can be written in the form 



dC 
~dt 
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Combining Eq. (PT|) and Eq. ()26|) and taking into account that the concen- 
tration of interstitial impurity atoms is essentially smaller than the concen- 
tration of substitutionally dissolved impurity atoms gives the final equation 
of impurity diffusion due to a kick-out mechanism 
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As can be seen from Eq.(|27p. the concentration of interstitial impurity 
atoms is not included in the explicit form in this equation, although Eq. (j27|) 
describes the diffusion due to a kick-out mechanism. The effective coefficients 
D{x) and v(x, x) of Eq.((27j) are the same as in Eq.(j2H), i.e. they repre- 
sent smooth and monotone functions of x- Moreover, the obtained equation 
Eq.(j2Zj) exactly matches the equation of impurity diffusion due to genera- 
tion, migration, and dissociation of the equilibrium pairs "impurity atom - 
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self-interstitial" , where the influence of elastic stress on the drift of the pairs 
is taken into account [Sni EZl • Thus, we can confirm the conclusion of [22] 
that macroscopic descriptions of impurity diffusion due to the kick-out and 
pair diffusion mechanisms are the same and analysis of the impurity concen- 
tration profiles does not allow to choose between these diffusion mechanisms. 
In our opinion, the difference can arise only in the case when the interstitial 
impurity atom makes one jump before conversion to the substitutional posi- 
tion. Because of this, it is possible to use a wide-spread software based on 
the models of pair diffusion for simulation of the diffusion processes due to a 
kick-out mechanism even if a stress-mediated migration occurs. 

5 Conclusions 

The diffusion equation for nonequilibrium interstitial impurity atoms, taking 
into account their different charge states and drift of all mobile interstitial 
species in the built-in electric field and in the field of elastic stresses, was 
obtained. The obtained generalized equation has the following characteristic 
features, i) This equation can describe the long-range migration of nonequi- 
librium impurity interstitials. ii) A set of diffusion equations written for 
the interstitial impurity atoms in each individual charge state is replaced 
by a single generalized equation, iii) The concentration of neutral impurity 
interstitials is the only variable that should be estimated in the obtained 
equation, iv) There is no explicit term describing the drift due to electric 
field in the equation obtained, v) Although the effective coefficients of the 
diffusion equation for interstitial impurity atoms represent nonlinear func- 
tions of the concentration of charged carriers, these functions are smooth 
and monotone. Due to these features, the derived equation is convenient for 
numerical solution. 

Based on the obtained generalized equation and mass action law for the 
interstitial impurity atoms being in equilibrium with substitutionally dis- 
solved impurity and self-interstitials, the equation of impurity diffusion due 
to a kick-out mechanism was derived. The obtained equation of impurity 
diffusion takes into account all charge states of interstitial impurity atoms 
and drift of interstitial species in the electric field and in the field of elastic 
stresses. The effective coefficients of this diffusion equation also represent 
smooth and monotone functions of the charge carrier concentration. More- 
over, the obtained equation exactly matches the equation of impurity diffu- 
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sion due to generation, migration, and dissociation of the equilibrium pairs 
"impurity atom - self-interstitial" , where the influence of elastic stresses on 
the drift of the pairs is taken into account. Thus, the conclusion of [22] 
concerning identity of the governing diffusion equations for both diffusion 
mechanisms is extended to the case of stress-mediated diffusion. 
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